Coulomb correlation effects in LaOFeAs: LDA+DMFT(QMC) study 
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Effects of Coulomb correlation on LaOFeAs electronic structure have been investigated by 
LDA+DMFT(QMC) method. The calculation results show that LaOFeAs is in the regime of in- 
termediate correlation strength with significant part of the spectral density moved from the Fermi 
energy to Hubbard bands. However the system is not on the edge of metal insulator-transition 
because increase of the Coulomb interaction parameter value from (7=4.0 eV to £7=5.0 eV did not 
result in insulator state. Correlations affect different d-orbitals not in the same way. ti g states 
(xz, yz and x 2 — y 2 orbitals) have higher energy due to crystal filed splitting and are nearly half- 
filled. Their spectral functions have pseudogap with Fermi energy position on the higher sub-band 
slope. Lower energy e s set of d-orbitals (3z 2 — r 2 and xy) have significantly larger occupancy values 
with typically metallic spectral functions. 

PACS numbers: 



Introduction. The very recent report of superconduc- 
tivity with the remarkable T c of 26K in LaOi-^F^FeAs 1 
has stimulated an intense experimental and theoretical 
activity. The first attempts to calculate electron-phonon 
coupling in this material^ have shown that it can not 
give so large superconducting transition temperature T c . 
That resembles the situation with high-T c cuprates and 
rise a question of the strength for Coulomb correlation 
effects in this material. Recently Haule et a$- have re- 
ported results of LDA+DMFT calculations for LaOFeAs 
with effective impurity model solved by continuous time 
quantum Monte-Carlo method 4 . Their conclusion is 
that this system is in strongly correlated regime on the 
edge of metal-insulator transition with only 20-30% of 
spectral weight in quasiparticle band. In the present 
work we report results of LDA+DMFT (QMC) calcula- 
tions for LaOFeAs with efective impurity model solved 
by standard quantum Monte-Carlo method in Hirsch- 
Fye algorithm 5 . Our results confirm that correlation 
effects in LaOFeAs indeed have sizable effect on elec- 
tronic structure. However the system is not on the edge 
of metal-insulator transition. While Haule et a£ claim 
that Coulomb parameter U increasing on only 0.5 eV 
from [7=4 eV value drives the system to insulating state, 
in our calculations even for [7=5.0 eV the solution is still 
clearly metallic. 

Computation details. Band structure of LaOFeAs was 
calculated in LDA 6 approximation by TB-LMTO-ASAl 
method. Crystal structure was taken from i . The results 
agree well with previously reported calculations by FP- 
LAPW method 8 . In Fig. Q] density of states for energy re- 
gion corresponding to the bands formed by Fe3d-orbitals 
is presented. Fermi level corresponding to 6 electrons per 
formula unit in d-band is situated near deep minimum 
between two sub-bands. These sub-bands have their ori- 
gin in bonding-antibonding separation due to strong hy- 



bridization between d-states in FeAs layers. 

Using eigenvalues and eigenfunctions from LDA calcu- 
lations for bands formed by Fe3d-orbitals we have cal- 
culated Wannier functions via projection procedur e 9 ' 10 . 
The calculated Hamiltonian matrix H(k) in Wannier 
functions basis has dimensions 10x10 (5 d-orbitals per 
Fe ion and two Fe ions in unit cell) and by construction 
has eigenvalues exactly the same as for ten LDA bands 
formed by Fe-3d orbitals. 

This Hamiltonian was used in LDA+DMFT— calcu- 
lations (for detailed description of the present computa- 
tion scheme seei -). The Coulomb interaction parame- 
ter value C/=4 eV and Hund's parameter J=0.7 eV used 
in our work were the same as in previous LDA+DMFT 
calculations by Haule and Kotliar— . The effective impu- 
rity model for DMFT was solved by QMC method in 
Hirsh-Fye algorithm 5 . Calculations were performed for 
the value of inverse temperature /3=10 eV^ 1 . Inverse 
temperature interval < r < j3 was divided in 100 slices. 
4 millions QMC sweeps were used in self-consistency loop 
within LDA+DMFT scheme and 12 millions of QMC 
sweeps were used to calculate spectral functions. 

Results and discussion. Total spectral function calcu- 
lated in LDA+DMFT in comparison with total density 
of states (DOS) obtained within LDA are presented in 
Fig. Q] (see also Fig. [3] for wider energy range). Correla- 
tion effects are strong enough to move part of the spectral 
density from Fermi energy to "shoulders" at -4 eV and 
+4eV. However in contrast to Haule et a£ spectra where 
strong Hubbard band at -4 eV is separated from the weak 
quasiparticle band by a deep wide depression in (-2 eV, 
eV) area, in our results a main part of spectral density 
is still in the interval from -2eV till +2 eV. In this en- 
ergy region LDA+DMFT spectral function preserves the 
general shape of LDA DOS with two subbands divided 
by a deep minimum. The essential result of correlation 
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FIG. 1: (Color online) DOS obtained within LDA approxi- 
mation (black line) and spectral function from LDA+DMFT 
for C/=4.0 eV, J=0.7 eV (color line). 



is the position of the Fermi energy that has moved from 
the edge of lower subband in LDA DOS to the top of 
the upper subband. Exactly on the Fermi energy a sharp 
peak has developed. 

LaOFeAs crystal structure has tetragonal symmetry 
and hence Green function for d-orbitals is diagonal with 
different matrix elements for every one of 3z 2 — r 2 , x 2 —y 2 , 
xy orbitals and doubly degenerate elements for xz, yz set 
of orbitals. Fe ion has tetrahedral coordination of four 
As ions with slight tetragonal distortion of the tetrahe- 
dron. In tetrahedral symmetry group Td five d-orbitals 
should be split by crystal field on low-energy doubly de- 
generate set 3z 2 — r 2 , xy corresponding to irreducible rep- 
resentation e^ 2 - and high-energy triply degenerate set 
x 2 — y 2 , xz, yz for representation t<z g . We have calculated 
Wannier functions energy and have found that t2 g -e g 
crystal field splitting parameter is very small A c f «0.25 
eV. The slight tetragonal distortion of the tetrahedron 
gives additional splitting of ti g and e g levels with the fol- 
lowing values for orbital energies (energy of the lowest 
xy orbital is taken as a zero): £^=0.0 eV, £3 Z 2_,.2=0.03 
eV, £^2,1,2=0.26 eV, e x 2_ y 2=QAl eV. Intersite hybridiza- 
tion between d-orbitals is much stronger than crystal field 
splitting (the total d-band width is s=s4 eV). As a result 
all five d-orbitals form common d-band with sub-bands 
originated not from crystal filed splitting, as it is the 
case in octahedral coordinated oxides, but due to strong 
bonding-antibonding separation (see Fig. [2]). rry-orbital 
that has lobes directed along Fe-Fe bond in FeAs planes 
shows specially strong bonding-antibonding separation. 

In Fig. [5] we present orbitally resolved spectral func- 
tions obtained in LDA+DMFT calculation with [7=4.0 
eV. Comparison with corresponding LDA partial DOS 
shows that different orbitals are affected by correlations 
not in the same way. x 2 — y 2 curve has nearly insula- 
tor shape with a pronounced pseudogap below the Fermi 
energy. Another t2 g orbitals xz, yz also show presence 
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FIG. 2: (Color online) PDOS obtained within LDA approxi- 
mation (green lines) and LDA+DMFT partial spectral func- 
tions for [7=4.0 eV, J=0.7 eV (blue lines). 



of peudogap but Fermi energy is on the slope of upper 
sub-band, with a general form typical to slightly electron 
doped Mott insulator. e g set of d-orbitals (3z 2 — r 2 and 
xy) show clearly metallic behaviour with well developed 
quasiparticle bands and no pseudogap. 

Occupation numbers obtained from LDA+DMFT 
calculations have the following values: 
yz,n zx =0.565, n 3z 2_ r 2=0.686 and n x ~ 



n 



-y- 



n xy =0.672, 

=0.512. These 



values mean that one of ti g states, x 2 —y 2 orbital that has 



highest crystal field energy (e x 2_ y 2—0Al eV) is nearly 
half-filled. Energy of the other of t 2g states (xz,yz- 
orbitals) is 0.15 eV lower (e xz ^ yz =Q.2& eV) and these or- 
bitals are slightly more occupied. e g set of d-orbitals 



(3z 2 



and xy) which have lowest energy due to crys- 



tal field splitting (s xy =0, £ 3z 2_ r 2=0.03 eV) have signifi- 
cantly larger occupancy. This fact agrees with our analy- 
sis of orbitally resolved spectral functions in Fig. [2j Half- 
filling of x 2 — y 2 orbital (n a .2_ y 2=0.512) corresponds to 
Mott insulator shape of the spectral function. Devia- 
tion from half-filling for xz,yz orbitals (n yz ,n zx — 0.565) 
results in light electron doped Mott insulator. The oc- 
cupancies of e g orbitals are much larger than t$ g values 
and corresponding spectral functions are characteristic 
for correlated metal. 

Authors of 3 report that "slightly enhanced Coulomb 
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FIG. 3: (Color online) Spectral function from LDA+DMFT 
for [7=4.0 eV (black line) and [7=5 eV (color line) ( J=0.7 eV 
for both). 



repulsion (U= 4.5 eV) opens the gap" so that the system 
is on the edge of metal-insulator transition. We have re- 
peated LDA+DMFT(QMC) calculations with even large 
Coulomb parameter [7=5.0 eV (Hund's parameter J=0.7 
eV) . The resulting total spectral function is presented in 
Fig. [3] in comparison with result for [7=4.0 eV. While 
increasing of U value leads to redistribution of spectral 
density from the peaks around Fermi to lower and upper 
Hubbard bands, the system still shows clearly metallic 
character with sharp peak on the Fermi energy. One can 
conclude that for [7=4.0 eV value the system is not close 
to metal-insulator transition. 

It is interesting to analyze evolution of orbitally re- 
solved spectral functions with increasing [7 value (Fig.HJ). 
For [7=5 eV t% g states comparing with [7=4 eV results 
(Fig. \2§ become fully insulating with quasiparticle peak 
absent for xz,yz orbitals and Fermi level position inside 
the gap for x 2 
d-orbitals (3z 2 — r 
tide peak on the Fermi energy. That type of electronic 
structure with partial localization resembles "Orbitally 
Selective Mott transition"— effect that was proposed for 



y spectral function. However e g set of 
2 and xy) still have a strong quasipar- 



Ca2- x Sr x Ru04 materials. Mott transition for x — y 
spectral function is not 100% pure because one can ob- 
serve weak peak on the Fermi energy inside the gap. This 
peak can be due to the mixing of x 2 — y 2 orbitals with 
metallic e g states. The evolution of orbital occupancies 
with increasing [7 value agrees with better localization 
of t 2g states for [7=5 eV: for (xz,yz),x 2 — y 2 orbitals 
occupancy values have decreased from (0.565, 0.512) to 
(0.535, 0.505) and for e g set of d-orbitals (3z 2 — r 2 and 
xy) the occupancies have increased from (0.686, 0.672) 
to (0.750, 0.672). 

We have found that inclusion of Hund parameter 
J=0.7 eV into calculations dramatically changes the ef- 
fect of correlations on LaOFeAs. LDA+DMFT(QMC) 
calculations with [7=4eV and J=0 eV results in very 
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FIG. 4: (Color online) PDOS obtained within LDA approxi- 
mation (green lines) and LDA+DMFT partial spectral func- 
tions for [7=5.0 eV, J=0.7 eV (blue lines). 



weakly correlated electronic structure with quasiparticle 
renormalization amplitude Z f«0.8. This agrees with the 
results o££, where effect of d-band degeneracy Nd on 
critical values for U/W ratio needed for metal-insulator 
transition was investigated. The authors oi^ have esti- 
mated the critical [7 value as U c « y/NgU® — NdJ, where 
U® is critical U value for non-degenerate Hubbard model. 
With LDA band width value W sw 4 eV and Coulomb in- 
teraction parameter [7= 4.0 eV their ratio U/W is close 
to 1. With J=0 orbital degeneracy N^—5 would lead to 
critical [7 value increased by a factor of y5 or effective 
decrease of U/W ratio more than twice thus resulting in 
weakly correlated regime. However J=0.7 eV will com- 
pensate increase of U c due to the y5 factor and U/W «1 
should result in intermediate correlation strength regime. 

In conclusion. The LDA-fDMFT(QMC) calcula- 
tions show that LaOFeAs is in intermediate correlation 
strength regime. Hubbard bands in spectral function are 
well pronounced but still a major part of spectral density 
is in the bands around Fermi energy. The system is not 
close to metal-insulator transition because even increase 
of Coulomb interaction parameter value on 1 eV does not 
lead to insulating state. 
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